load('sandiantu.mat')
x = zeros(1,150);
y = zeros(1,150);
z = zeros(1,150);
temt = 1;
for q = 1:10
    for w = 1:5
        for r = 1:5
            x(temt) = q;
            y(temt) = w;
            z(temt) = all_rate(q,w,r);
            temt = temt +1 ;
        end
    end
end
scatter3(x,y,z,40,z,'filled')
ax = gca ;
ax.XDir = 'reverse';
view(-31,14)
xlabel('cap')
ylabel('r0')
zlabel('Rick_w sss')




        